Synergy of Small Antiviral Molecules on a Black-Phosphorus Nanocarrier: Machine Learning and Quantum Chemical Simulation Insights.

Favipiravir (FP) and Ebselen (EB) belong to a broad range of antiviral drugs that have shown active potential as medications against many viruses. Employing molecular dynamics simulations and machine learning (ML) combined with van der Waals density functional theory, we have uncovered the binding characteristics of these two antiviral drugs on a phosphorene nanocarrier. Herein, by using four different machine learning models (i.e., Bagged Trees, Gaussian Process Regression (GPR), Support Vector Regression (SVR), and Regression Trees (RT)), the Hamiltonian and the interaction energy of antiviral molecules in a phosphorene monolayer are trained in an appropriate way. However, training efficient and accurate models for approximating the density functional theory (DFT) is the final step in using ML to aid in the design of new drugs. To improve the prediction accuracy, the Bayesian optimization approach has been employed to optimize the GPR, SVR, RT, and BT models. Results revealed that the GPR model obtained superior prediction performance with an R2 of 0.9649, indicating that it can explain 96.49% of the data's variability. Then, by means of DFT calculations, we examine the interaction characteristics and thermodynamic properties in a vacuum and a continuum solvent interface. These results illustrate that the hybrid drug is an enabled, functionalized 2D complex with vigorous thermostability. The change in Gibbs free energy at different surface charges and temperatures implies that the FP and EB molecules are allowed to adsorb from the gas phase onto the 2D monolayer at different pH conditions and high temperatures. The results reveal a valuable antiviral drug therapy loaded by 2D biomaterials that may possibly open a new way of auto-treating different diseases, such as SARS-CoV, in primary terms.

Insighting the optoelectronic, charge transfer and biological potential of benzo-thiadiazole and its derivatives.

The current investigation applies the dual approach containing quantum chemical and molecular docking techniques to explore the potential of benzothiadiazole (BTz) and its derivatives as efficient electronic and bioactive materials. The charge transport, electronic and optical properties of BTz derivatives are explored by quantum chemical techniques. The density functional theory (DFT) and time dependent DFT (TD-DFT) at B3LYP/6-31G** level of theory utilized to optimize BTz and newly designed ligands at the ground and first excited states, respectively. The heteroatoms substitution effects on different properties of 4,7-bis(4-methylthiophene-2yl) benzo[c] [1,2,5]thiadiazole (BTz2T) as initial compound are studied at molecular level. Additionally, we also study the possible inhibition potential of COVID-19 from benzothiadiazole (BTz) containing derivatives by implementing the grid based molecular docking methods. All the newly designed ligands docked with the main protease (MPRO:PDB ID 6LU7) protein of COVID-19 through molecular docking methods. The studied compounds showed strong binding affinities with the binding site of MPRO ranging from -6.9 to -7.4 kcal/mol. Furthermore, the pharmacokinetic properties of the ligands are also studied. The analysis of these results indicates that the studied ligands might be promising drug candidates as well as suitable for photovoltaic applications.

Fabrication, DFT Calculation, and Molecular Docking of Two Fe(III) Imine Chelates as Anti-COVID-19 and Pharmaceutical Drug Candidate.

Two tetradentate dibasic chelating Schiff base iron (III) chelates were prepared from the reaction of 2,2'-((1E,1'E)-(1,2-phenylenebis(azanylylidene))bis(methanylylidene))bis(4-bromophenol) (PDBS) and 2,2'-((1E,1'E)-((4-chloro-1,2-phenylene)bis(azanylylidene))-bis(methanylylidene))bis(4-bromophenol) (CPBS) with Fe3+ ions. The prepared complexes were fully characterized with spectral and physicochemical tools such as IR, NMR, CHN analysis, TGA, UV-visible spectra, and magnetic moment measurements. Moreover, geometry optimizations for the synthesized ligands and complexes were conducted using the Gaussian09 program through the DFT approach, to find the best structures and key parameters. The prepared compounds were tested as antimicrobial agents against selected strains of bacteria and fungi. The results suggests that the CPBSFe complex has the highest activity, which is close to the reference. An MTT assay was used to screen the newly synthesized compounds against a variety of cell lines, including colon cancer cells, hepatic cellular carcinoma cells, and breast carcinoma cells. The results are expressed by IC50 value, in which the 48 µg/mL value of the CPBSFe complex indicates its success as a potential anticancer agent. The antioxidant behavior of the two imine chelates was studied by DPPH assay. All the tested imine complexes show potent antioxidant activity compared to the standard Vitamin C. Furthermore, the in vitro assay and the mechanism of binding and interaction efficiency of the tested samples with the receptor of COVID-19 core protease viral protein (PDB ID: 6lu7) and the receptor of Gram-negative bacteria (Escherichia coli, PDB ID: 1fj4) were investigated using molecular docking experiments.

A Competition between Hydrogen, Stacking, and Halogen Bonding in N-(4-((3-Methyl-1,4-dioxo-1,4-dihydronaphthalen-2-yl)selanyl)phenyl)acetamide: Structure, Hirshfeld Surface Analysis, 3D Energy Framework Approach, and DFT Calculation.

N-(4-((3-Methyl-1,4-dioxo-1,4-dihydronaphthalen-2-yl)selanyl)phenyl)acetamide (5), C19H15NO3Se, was prepared in two steps from 4,4'-diselanediyldianiline (3) via reduction and subsequent nucleophilic reaction with 2-methyl-3-bromo-1,4-naphthalenedione, followed by acetylation with acetic anhydride. The cytotoxicity was estimated against 158N and 158JP oligodendrocytes and the redox profile was also evaluated using different in vitro assays. The technique of single-crystal X-ray diffraction is used to confirm the structure of compound 5. The enantiopure 5 crystallizes in space group P21 with Flack parameter 0.017 (8), exhibiting a chiral layered absolute structure. Molecular structural studies showed that the crystal structure is foremost stabilized by N-H···O and relatively weak C-H···O contacts between molecules, and additionally stabilized by weak C-H···π and Se···N interactions. Hirshfeld surface analysis is used to quantitatively investigate the noncovalent interactions that stabilize crystal packing. Framework energy diagrams were used to graphically represent the stabilizing interaction energies for crystal packing. The analysis of the energy framework shows that the interactions energies of and C-H···π and C-O···π are primarily dispersive and are the crystal's main important forces. Density functional theory (DFT) calculations were used to determine the compound's stability, chemical reactivity, and other parameters by determining the HOMO-LUMO energy differences. The determination of its optimized surface of the molecular electrostatic potential (MEP) was also carried out. This study was conducted to demonstrate both the electron-rich and electron-poor sites.

Jusanin, a New Flavonoid from Artemisia commutata with an In Silico Inhibitory Potential against the SARS-CoV-2 Main Protease.

A new flavonoid, Jusanin, (1) has been isolated from the aerial parts of Artemisia commutata. The chemical structure of Jusanin has been elucidated using 1D, 2D NMR, and HR-Ms spectroscopic methods to be 5,2',4'-trihydroxy-6,7,5'-trimethoxyflavone. Being new in nature, the inhibition potential of 1 has been estimated against SARS-CoV-2 using different in silico techniques. Firstly, molecular similarity and fingerprint studies have been conducted for Jusanin against co-crystallized ligands of eight different SARS-CoV-2 essential proteins. The studies indicated the similarity between 1 and X77, the co-crystallized ligand SARS-CoV-2 main protease (PDB ID: 6W63). To confirm the obtained results, a DFT study was carried out and indicated the similarity of (total energy, HOMO, LUMO, gap energy, and dipole moment) between 1 and X77. Accordingly, molecular docking studies of 1 against the target enzyme have been achieved and showed that 1 bonded correctly in the protein's active site with a binding energy of -19.54 Kcal/mol. Additionally, in silico ADMET in addition to the toxicity evaluation of Jusanin against seven models have been preceded and indicated the general safety and the likeness of Jusanin to be a drug. Finally, molecular dynamics simulation studies were applied to investigate the dynamic behavior of the Mpro-Jusanin complex and confirmed the correct binding at 100 ns. In addition to 1, three other metabolites have been isolated and identified to be сapillartemisin A (2), methyl-3-[S-hydroxyprenyl]-cumarate (3), and ß-sitosterol (4).

Ab Initio Insight into the Interaction of Metal-Decorated Fluorinated Carbon Fullerenes with Anti-COVID Drugs.

We theoretically investigated the adsorption of two common anti-COVID drugs, favipiravir and chloroquine, on fluorinated C60 fullerene, decorated with metal ions Cr3+, Fe2+, Fe3+, Ni2+. We focused on the effect of fluoridation on the interaction of fullerene with metal ions and drugs in an aqueous solution. We considered three model systems, C60, C60F2 and C60F48, and represented pristine, low-fluorinated and high-fluorinated fullerenes, respectively. Adsorption energies, deformation of fullerene and drug molecules, frontier molecular orbitals and vibrational spectra were investigated in detail. We found that different drugs and different ions interacted differently with fluorinated fullerenes. Cr3+ and Fe2+ ions lead to the defluorination of low-fluorinated fullerenes. Favipiravir also leads to their defluorination with the formation of HF molecules. Therefore, fluorinated fullerenes are not suitable for the delivery of favipiravir and similar drugs molecules. In contrast, we found that fluorine enhances the adsorption of Ni2+ and Fe3+ ions on fullerene and their activity to chloroquine. Ni2+-decorated fluorinated fullerenes were found to be stable and suitable carriers for the loading of chloroquine. Clear shifts of infrared, ultraviolet and visible spectra can provide control over the loading of chloroquine on Ni2+-doped fluorinated fullerenes.

Metal oxide nanocage as drug delivery systems for Favipiravir, as an effective drug for the treatment of COVID-19: a computational study.

This paper is a summary of research that looks at the potential of fullerene-like (MO)12 nanoclusters (NCs) in drug-carrying systems using density functional theory. Favipiravir/Zn12O12 (- 34.80 kcal/mol), Favipiravir/Mg12O12 (- 34.98 kcal/mol), and Favipiravir/Be12O12 (- 30.22 kcal/mol) were rated in order of drug adsorption degrees. As a result, Favipiravir attachment to (MgO)12 and (ZnO)12 might be simple, increasing Favipiravir loading efficiency. In addition, the quantum theory of atoms in molecules (QTAIM) assessment was utilized to look at the interactions between molecules. The FMO, ESP, NBO, and Eads reactivity patterns were shown to be in excellent agreement with the QTAIM data. The electrostatic properties of the system with the biggest positive charge on the M atom and the largest Eads were shown to be the best. This system was shown to be the best attraction site for nucleophilic agents. The findings show that (MgO)12 and (ZnO)12 have great carrier potential and may be used in medication delivery.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex.

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

Applications of density functional theory in COVID-19 drug modeling.

The rapidly evolving Coronavirus 2019 (COVID-19) pandemic has led to millions of deaths around the world, highlighting the pressing need to develop effective antiviral pharmaceuticals. Recent efforts with computer-aided rational drug discovery have allowed detailed examination of drug-macromolecule interactions primarily by molecular mechanics (MM) techniques. Less widely applied in COVID-19 drug modeling is density functional theory (DFT), a quantum mechanics (QM) method that enables electronic structure calculations and elucidations of reaction mechanisms. Here, we review recent advances in applying DFT in molecular modeling studies of COVID-19 pharmaceuticals. We start by providing an overview of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) drugs and targets, followed by a brief introduction to DFT. We then provide a discussion of different approaches by which DFT has been applied. Finally, we discuss essential factors to consider when incorporating DFT in future drug modeling research.

DFT calculations to investigate silver ions as a virucide from SARS-CoV-2.

The world has face the COVID-19 pandemic which has already caused millions of death. Due to the urgency in fighting the virus, we study five residues of free amino acids present in the structure of the SARS-CoV-2 spike protein (S). We investigated the spontaneous interaction between amino acids and silver ions (Ag+), considering these ions as a virucide chemical agent for SARS-CoV-2. The amino acid-Ag+ systems were investigated in a gaseous medium and a simulated water environment was described with a continuum model (PCM) the calculations were performed within the framework of density functional theory (DFT). Calculations related to the occupied orbitals of higher energy showed that Ag+ has a tendency to interact with the nitrile groups (-NH). The negative values of the Gibbs free energies show that the interaction process between amino acids-Ag+ in both media occurs spontaneously. There is a decrease in Gibbs free energy from the amino acid-Ag+ interactions immersed in a water solvation simulator.

In Silico Exploration of Potential Natural Inhibitors against SARS-Cov-2 nsp10.

In continuation of our previous effort, different in silico selection methods were applied to 310 naturally isolated metabolites that exhibited antiviral potentialities before. The applied selection methods aimed to pick the most relevant inhibitor of SARS-CoV-2 nsp10. At first, a structural similarity study against the co-crystallized ligand, S-Adenosyl Methionine (SAM), of SARS-CoV-2 nonstructural protein (nsp10) (PDB ID: 6W4H) was carried out. The similarity analysis culled 30 candidates. Secondly, a fingerprint study against SAM preferred compounds 44, 48, 85, 102, 105, 182, 220, 221, 282, 284, 285, 301, and 302. The docking studies picked 48, 182, 220, 221, and 284. While the ADMET analysis expected the likeness of the five candidates to be drugs, the toxicity study preferred compounds 48 and 182. Finally, a density-functional theory (DFT) study suggested vidarabine (182) to be the most relevant SARS-Cov-2 nsp10 inhibitor.

Synthesis, In Silico Prediction and In Vitro Evaluation of Antimicrobial Activity, DFT Calculation and Theoretical Investigation of Novel Xanthines and Uracil Containing Imidazolone Derivatives.

Novel xanthine and imidazolone derivatives were synthesized based on oxazolone derivatives 2a-c as a key intermediate. The corresponding xanthine 3-5 and imidazolone derivatives 6-13 were obtained via reaction of oxazolone derivative 2a-c with 5,6-diaminouracils 1a-e under various conditions. Xanthine compounds 3-5 were obtained by cyclocondensation of 5,6-diaminouracils 1a-c with different oxazolones in glacial acetic acid. Moreover, 5,6-diaminouracils 1a-e were reacted with oxazolones 2a-c in presence of drops of acetic acid under fused condition yielding the imidazolone derivatives 6-13. Furthermore, Schiff base of compounds 14-16 were obtained by condensing 5,6-diaminouracils 1a,b,e with 4-dimethylaminobenzaldehyde in acetic acid. The structural identity of the resulting compounds was resolved by IR, 1H-, 13C-NMR and Mass spectral analyses. The novel synthesized compounds were screened for their antifungal and antibacterial activities. Compounds 3, 6, 13 and 16 displayed the highest activity against Escherichia coli as revealed from the IC50 values (1.8-1.9 µg/mL). The compound 16 displayed a significant antifungal activity against Candia albicans (0.82 µg/mL), Aspergillus flavus (1.2 µg/mL) comparing to authentic antibiotics. From the TEM microgram, the compounds 3, 12, 13 and 16 exhibited a strong deformation to the cellular entities, by interfering with the cell membrane components, causing cytosol leakage, cellular shrinkage and irregularity to the cell shape. In addition, docking study for the most promising antimicrobial tested compounds depicted high binding affinity against acyl carrier protein domain from a fungal type I polyketide synthase (ACP), and Baumannii penicillin- binding protein (PBP). Moreover, compound 12 showed high drug- likeness, and excellent pharmacokinetics, which needs to be in focus for further antimicrobial drug development. The most promising antimicrobial compounds underwent theoretical investigation using DFT calculation.

Advancing insights on ß-cyclodextrin inclusion complexes with SSRIs through lens of X-ray diffraction and DFT calculation.

Depression-the global crisis hastened by the coronavirus outbreak, can be efficaciously treated by the selective serotonin reuptake inhibitors (SSRIs). Cyclodextrin (CD) inclusion complexation is a method of choice for reducing side effects and improving bioavailability of drugs. Here, we investigate in-depth the ß-CD encapsulation of sertraline (STL) HCl (1) and fluoxetine (FXT) HCl (2) by single-crystal X-ray diffraction and DFT complete-geometry optimization, in comparison to the reported complex of paroxetine (PXT) base. X-ray analysis unveiled the 2:2 ß-CD-STL/FXT complexes with two drug molecules inserting their halogen-containing aromatic ring in the ß-CD dimeric cavity, which are stabilized by the interplay of intermolecular O2-H⋯N1-H⋯O3 H-bonds, C3/C5-H⋯π and halogen⋯halogen interactions. Similarly, the 1:1 ß-CD-tricyclic-antidepressant (TCA) complexes have an exclusive inclusion mode of the aromatic ring, which is maintained by C3/C5-H⋯π interactions. By contrast, the 2:1 ß-CD-PXT complex has a total inclusion that is stabilized by host-guest O6-H⋯N1-H⋯O5 H-bonds and C3-H⋯π interactions. The inherent stabilization energies of 1 and 2 evaluated using DFT calculation suggested that the improved thermodynamic stabilities via CD encapsulation facilitates the reduction of drug side effects. Moreover, the SSRI conformational flexibilities are thoroughly discussed for understanding of their pharmacoactivity.

Quantum simulations of SARS-CoV-2 main protease Mpro enable high-quality scoring of diverse ligands.

The COVID-19 pandemic has led to unprecedented efforts to identify drugs that can reduce its associated morbidity/mortality rate. Computational chemistry approaches hold the potential for triaging potential candidates far more quickly than their experimental counterparts. These methods have been widely used to search for small molecules that can inhibit critical proteins involved in the SARS-CoV-2 replication cycle. An important target is the SARS-CoV-2 main protease Mpro, an enzyme that cleaves the viral polyproteins into individual proteins required for viral replication and transcription. Unfortunately, standard computational screening methods face difficulties in ranking diverse ligands to a receptor due to disparate ligand scaffolds and varying charge states. Here, we describe full density functional quantum mechanical (DFT) simulations of Mpro in complex with various ligands to obtain absolute ligand binding energies. Our calculations are enabled by a new cloud-native parallel DFT implementation running on computational resources from Amazon Web Services (AWS). The results we obtain are promising: the approach is quite capable of scoring a very diverse set of existing drug compounds for their affinities to M pro and suggest the DFT approach is potentially more broadly applicable to repurpose screening against this target. In addition, each DFT simulation required only ~ 1 h (wall clock time) per ligand. The fast turnaround time raises the practical possibility of a broad application of large-scale quantum mechanics in the drug discovery pipeline at stages where ligand diversity is essential.

Overestimation of 3α- over 3ß-25-Hydroxyvitamin D3 Levels in Serum: A Mechanistic Rationale for the Different Mass Spectral Properties of the Vitamin D Epimers.

The metabolism of vitamin D3 includes a parallel C-3 epimerization pathway-in addition to the standard metabolic processes for vitamin D3-reversing the stereochemical configuration of the -OH group at carbon-3 (ß→α). While the biological function of the 3α epimer has not been elucidated yet, the additional species cannot be neglected in the analytical determination of vitamin D3, as it has the potential to introduce analytical errors if not properly accounted for. Recently, some inconsistent mass spectral behavior was seen for the 25-hydroxyvitamin D3 (25(OH)D3) epimers during quantification using electrospray LC-MS/MS. The present work extends that of Flynn et al. ( Ann. Clin. Biochem. 2014, 51, 352-559) and van den Ouweland et al. ( J. Chromatogr. B 2014, 967, 195-202), who reported larger electrospray ionization response factors for the 3α epimer of 25(OH)D3 in human serum samples as compared to the regular 3ß variant. The present work was concerned with the mechanistic reasons for these differences. We used a combination of electrospray ionization, atmospheric pressure chemical ionization, and density functional theory calculations to uncover structural dissimilarities between the epimers. A plausible mechanism is described based on intramolecular hydrogen bonding in the gas phase, which creates a small difference of proton affinities between the epimers. More importantly, this mechanism allows the explanation of the different ionization efficiencies of the epimers based on kinetic control of the ionization process, where ionization initially takes place at the hydroxyl group with subsequent proton transfer to a basic carbon atom. The barrier for this transfer differs between the epimers and is in direct competition with H2O elimination from the protonated hydroxyl group. The "hidden" site of high gas phase basicity was revealed through computational calculations and appears to be inaccessible via direct protonation.

An Androsterone-H2 @C60 hybrid: Synthesis, Properties and Molecular Docking Simulations with SARS-Cov-2.

We report the synthesis and characterization of a fullerene-steroid hybrid that contains H2 @C60 and a dehydroepiandrosterone moiety synthesized by a cyclopropanation reaction with 76 % yield. Theoretical calculations at the DFT-D3(BJ)/PBE 6-311G(d,p) level predict the most stable conformation and that the saturation of a double bond is the main factor causing the upfield shielding of the signal appearing at -3.13 ppm, which corresponds to the H2 located inside the fullerene cage. Relevant stereoelectronic parameters were also investigated and reinforce the idea that electronic interactions must be considered to develop studies on chemical-biological interactions. A molecular docking simulation predicted that the binding energy values for the protease-hybrid complexes were -9.9 kcal/mol and -13.5 kcal/mol for PLpro and 3CLpro respectively, indicating the potential use of the synthesized steroid-H2 @C60 as anti-SARS-Cov-2 agent.

Pt-, Rh-, Ru-, and Cu-Single-Wall Carbon Nanotubes Are Exceptional Candidates for Design of Anti-Viral Surfaces: A Theoretical Study.

As SARS-CoV-2 is spreading rapidly around the globe, adopting proper actions for confronting and protecting against this virus is an essential and unmet task. Reactive oxygen species (ROS) promoting molecules such as peroxides are detrimental to many viruses, including coronaviruses. In this paper, metal decorated single-wall carbon nanotubes (SWCNTs) were evaluated for hydrogen peroxide (H2O2) adsorption for potential use for designing viral inactivation surfaces. We employed first-principles methods based on the density functional theory (DFT) to investigate the capture of an individual H2O2 molecule on pristine and metal (Pt, Pd, Ni, Cu, Rh, or Ru) decorated SWCNTs. Although the single H2O2 molecule is weakly physisorbed on pristine SWCNT, a significant improvement on its adsorption energy was found by utilizing metal functionalized SWCNT as the adsorbent. It was revealed that Rh-SWCNT and Ru-SWCNT systems demonstrate outstanding performance for H2O2 adsorption. Furthermore, we discovered through calculations that Pt- and Cu-decorated SWNCT-H2O2 systems show high potential for filters for virus removal and inactivation with a very long shelf-life (2.2 × 1012 and 1.9 × 108 years, respectively). The strong adsorption of metal decorated SWCNTs and the long shelf-life of these nanomaterials suggest they are exceptional candidates for designing personal protection equipment against viruses.

Key Interacting Residues between RBD of SARS-CoV-2 and ACE2 Receptor: Combination of Molecular Dynamics Simulation and Density Functional Calculation.

The spike protein of SARS-CoV-2 binds to the ACE2 receptor via its receptor-binding domain (RBD), with the RBD-ACE2 complex presenting an essential molecular target for vaccine development to stall the virus infection proliferation. The computational analyses at molecular, amino acid (AA), and atomic levels have been performed systematically to identify the key interacting AAs in the formation of the RBD-ACE2 complex for SARS-CoV and SARS-CoV-2 with its Alpha and Beta variants. Our study uses the molecular dynamics (MD) simulations with the molecular mechanics generalized Born surface area (MM-GBSA) method to predict the binding free energy (BFE) and to determine the actual interacting AAs, as well as two ab initio quantum chemical protocols based on the density functional theory (DFT) implementation. Based on MD results, Q493, Y505, Q498, N501, T500, N487, Y449, F486, K417, Y489, F456, Y495, and L455 have been identified as hotspots in SARS-CoV-2 RBD, while those in ACE2 are K353, K31, D30, D355, H34, D38, Q24, T27, Y83, Y41, and E35. RBD with Alpha and Beta variants has slightly different interacting AAs due to N501Y mutation. Both the electrostatic and hydrophobic interactions are the main driving force to form the AA-AA binding pairs. We confirm that Q493, Q498, N501, F486, K417, and F456 in RBD are the key residues responsible for the tight binding of SARS-CoV-2 with ACE2 compared to SARS-CoV. RBD with the Alpha variant binds with ACE2 stronger than the wild-type RBD or Beta. In the Beta variant, K417N reduces the binding, E484K slightly enhances it, and N501Y significantly increases it as in Alpha. The DFT results reveal that N487, Q493, Y449, T500, G496, G446, and G502 in RBD of SARS2 form pairs via specific hydrogen bonding with Q24, H34, E35, D38, Y41, Q42, and K353 in ACE2.

Virtual Screening of Phytochemicals by Targeting HR1 Domain of SARS-CoV-2 S Protein: Molecular Docking, Molecular Dynamics Simulations, and DFT Studies.

The recent COVID-19 pandemic has impacted nearly the whole world due to its high morbidity and mortality rate. Thus, scientists around the globe are working to find potent drugs and designing an effective vaccine against COVID-19. Phytochemicals from medicinal plants are known to have a long history for the treatment of various pathogens and infections; thus, keeping this in mind, this study was performed to explore the potential of different phytochemicals as candidate inhibitors of the HR1 domain in SARS-CoV-2 spike protein by using computer-aided drug discovery methods. Initially, the pharmacological assessment was performed to study the drug-likeness properties of the phytochemicals for their safe human administration. Suitable compounds were subjected to molecular docking to screen strongly binding phytochemicals with HR1 while the stability of ligand binding was analyzed using molecular dynamics simulations. Quantum computation-based density functional theory (DFT) analysis was constituted to analyze the reactivity of these compounds with the receptor. Through analysis, 108 phytochemicals passed the pharmacological assessment and upon docking of these 108 phytochemicals, 36 were screened passing a threshold of -8.5 kcal/mol. After analyzing stability and reactivity, 5 phytochemicals, i.e., SilybinC, Isopomiferin, Lycopene, SilydianinB, and Silydianin are identified as novel and potent candidates for the inhibition of HR1 domain in SARS-CoV-2 spike protein. Based on these results, it is concluded that these compounds can play an important role in the design and development of a drug against COVID-19, after an exhaustive in vitro and in vivo examination of these compounds, in future.

A Combined Molecular Docking and Density Functional Theory Nuclear Magnetic Resonance Study of Trans-Dehydrocrotonin Interacting with COVID-19 Main Protease and Severe Acute Respiratory Syndrome Coronavirus 2 3C-Like Protease.

For the development of drugs that treat SARS-CoV-2, the fastest way is to find potential molecules from drugs already on the market. Unfortunately, there is currently no specific drug or treatment for COVID-19. Among all structural proteins in SARS-CoV, the spike protein is the main antigenic component responsible for inducing host immune responses, neutralizing antibodies, and/or protecting immunity against virus infection. Molecular docking is a technique used to predict whether a molecule will bind to another. It is usually a protein to another or a protein to a binding compound. Natural products are potential binders in several studies involving coronavirus. The structure of the ligand plays a fundamental role in its biological properties. The nuclear magnetic resonance technique is one of the most powerful tools for the structural determination of ligands from the origin of natural products. Nowadays, molecular modeling is an important accessory tool to experimentally got nuclear magnetic resonance data. In the present work, molecular docking studies aimed is to investigate the limiting affinities of trans-dehydrocrotonin molecule and to identify the main amino acid residues that could play a fundamental role in their mechanism of action of the SARS-CoV spike protein. Another aim of this work is all about to evaluate 10 hybrid functionalities, along with three base pairs using computational programs to discover which ones are more reliable with the experimental result the best computational method to study organic compounds. We compared the results between the mean absolute deviation (MAD) and root-mean-square deviation (RMSD) of the molecules, and the smallest number between them was the best result. The positions assumed by the ligands in the active site of the spike glycoprotein allow assuming associations with different local amino acids.